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Abstract. We present numerical simulations of the SU(2) Boltzmann-Vlasov 
^ i' equation including both hard elastic particle collisions and soft interactions mediated 



by classical Yang-Mills fields. We provide an estimate of the coupling of jets to a hot 

(J_) ' isotropic plasma which is independent of infrared cutoffs. In addition, we investigate jet 

j-^^' propagation in anisotropic plasmas, as created in heavy-ion collisions. The broadening 

of jets is found to be stronger along the beam line than in azimuth, due to the creation 

of field configurations with B± > E± and Ez > Bz via plasma instabilities. 
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High transverse niomentum jets produced in heavy-ion collisions represent a valuable 
tool for studies of the properties of the hot parton plasma produced in the central 
rapidity region [T]. However, present estimates of the strength of the coupling of jets to 
a QCD plasma are sensitive to infrared cutoffs. We employ a numerical simulation of 
the Boltzmann-Vlasov equation, which is coupled to the Yang-Mills equation for the soft 
gluon degrees of freedom. Soft momentum exchanges between particles are mediated by 
the fields, while hard momentum exchanges are described by a collision term including 
binary elastic collisions. This way, we are able to provide an estimate of the coupling of 
jets to a hot plasma which is independent of infrared cutoffs. 

The longitudinal expansion of the plasma may lead to a strongly anisotropic 
momentum distribution in the local rest frame, during the very early stages of the 
plasma evolution. Due to anisotropies in the particle momentum distributions plasma 
instabilities appear [2]. These lead to the formation of long- wavelength chromo-fields 
with E^ > Bz and B± > E±, which affect the propagation of a hard jet and of its induced 
hard radiation field. This may provide an explanation for the observed asymmetry in 
measurements of dihadron-correlations. Here a much stronger broadening of jets in 
pseudorapidity (rj) than in azimuthal angle (0) has been observed [H [3] . 
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2. Boltzmann-Vlasov equation for non-Abelian gauge theories 

We solve the classical transport equation for hard gluons with SU(2) color charge 
including hard binary collisions 

t/^ [d, + gq^'F^^d; + gr'^A'^{x)q^d,.] f = C, (1) 

where / = f{x,p,q) denotes the single-particle phase space distribution. It is coupled 
self-consistently to the Yang-Mills equation for the soft gluon fields. The collision term 
contains all binary collisions, descibed by the leading-order gg —>■ gg tree-level diagrams. 
We replace the distribution f{x,p,q) by a large number of test particles, which 
leads to Wong's equations [1] 

Mt) = v,(t) , Mt) = gqtit) {E%t) + v,(t) x B'^(t)) , 

g,(t) =~igvt{t)[A^{t),q,{t)], (2) 

for the i-th test particle, whose coordinates are Xj(t), Pi{t), and qf(t). The time evolution 
of the Yang-Mills field is determined by the standard Hamiltonian method [5] in A^ = 
gauge. See [6l EJ [HI [9] for more details. The collision term is incorporated using the 
stochastic method [lO]. The total cross section is given by a2^2 = J^L ^dq^ ■, where 
we have introduced a lower cutoff k* . To avoid double-counting, this cutoff should be 
on the order of the hardest field mode that can be represented on the given lattice, 
k* ~ vr/a, with the lattice spacing a. For a more detailed discussion see [9]. 

3. Jet broadening in an isotropic plasma 

We first consider a heat-bath of particles with a density of Ug = 10/fm^ and an average 
particle momentum of 3T = 12 GeV. For a given lattice (resp. k*) we take the initial 
energy density of the thermalized fields to be j dP'k / {2-n)^ kf-BQsc{k)Q{k* — k), where 
fBoseik) = ng/{2T^((3))/{e''^'^ — 1) is a Bose distribution normalized to the assumed 
particle density n^, and ( is the Riemann zeta function. The initial spectrum is fixed 
to Coulomb gauge and Ai ~ 1/k. We measure the momentum broadening {p\){t) 
of high-energy test particles ip/^T ^ 5) passing through this medium. Fig. [T] shows 
that in the coUisionless case, C = 0, the broadening is stronger on larger lattices 
which accommodate harder field modes. However, Fig. [2] demonstrates that collisions 
with momentum exchange larger than fc*(a) compensate for this growth and lead to 
approximately lattice-spacing independent results. 

A related transport coefficient is q [llj. It is the typical momentum transfer 
(squared) per collision divided by the mean-free path, which is nothing but {p\)(t)/t. 
From Fig. El g ^ 2.2 GeVVfm for N^ = 2, n^ = lO/fm^ and p/(3T) ^ 5. We have 
verified that q does not depend on the temperature T as long as the particle density Ug 
and the ratio p/T is fixed. Due to the independence of q of the temperature and its 
proportionality to the density n, we can scale to physical densities for a QGP created 
at RHIC. We adjust for the different color factors in SU(3), and find q ~ 5.6 GeV^/fm, 
at T = 400 MeV, E^et ~ 20 GeV {p/ST = 16) in a system of quarks and gluons. 
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Figure 1. Momentum difFusion 
caused by particle-field interac- 
tions only. 



Figure 2. Momentum difFusion 
by both particle-field and direct 
particle-particle interactions. 
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Figure 3. Momentum broaden- 
ing of a jet transverse to its ini- 
tial momentum. 



Figure 4. Ratios of field energy 
densities. 



4. Jet broadening in an unstable plasma 

In heavy-ion collisions, locally anisotropic momentum distributions may emerge due to 
the longitudinal expansion. Such anisotropies generically give rise to instabilities [SllTlE]- 
Here, we investigate their effect on the momentum broadening of jets, including the ef- 
fect of collisions. The initial anisotropic momentum distribution for the hard plasma 

gluons is taken to be /(p) = ng{ — \ (^(px) exp(— p_|_/ph) , with p^ = yfpl + p^- We 
initialize small-amplitude fields sampled from a Gaussian distribution and set k* ~ ph, 
for the reasons alluded to above. We add additional high momentum particles with 
Px = 12 ph and px = 6ph, respectively, to investigate the broadening in the y and z di- 
rections via the variances K±{px) '■= ^((Aj9_|_)^) , Kz{px) '■= -^{{^PzY) ■ The ratio Kz/k± 
can be roughly associated with the ratio of jet correlation widths in azimuth and ra- 
pidity: K,z/i^± ~ {^v)/{^4')- Experimental data on dihadron correlation functions for 
central Au-|-Au collisions at y/s = 200 GeV p] are consistent with k^/kj, ~ 3 [T2]. Fig. [3] 
shows the time evolution of (p^) and of (p^). During the period of instability and for 
both jet energies we find Hz/i^± ~ 2.3 . The explanation for the larger broadening along 
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the beam axis is as follows. In the Abelian case the instability generates predominantly 
transverse magnetic fields which defiect the particles in the z-direction [T3]. Although 
the interactions are a lot less trivial, in a non-Abelian plasma the instability creates 
large domains of strong chromo-electric and -magnetic fields with E^ > B^, aside from 
B± > E±_ (Fig.Hj). The field configurations are such that particles are defiected preferen- 
tially in the longitudinal z-direction (to restore isotropy). Fig. [5]shows the filamentation 
of the current and the domains of magnetic fields generated by the instability. 



t= 12.4826 m, 




Figure 5. Slices in the x-z-plane at fixed y = L/2 of the current in the x-direction, 
Jx, and the three color components of the chromo-magnetic field in the j/-direction. 
Filaments are nicely visible. Scales in lattice units. to 5 • 10^*^ for the current, 
—4 • 10~^ to 4 • 10~^ for the chromo-magnetic fields. 
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